Accretion Disk Instabilities, CDM models and their role in 

Quasar Evolution 
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ABSTRACT 

We have developed a consistent analytical model to describe the observed evolution of the 
quasar luminosity function. Our model combines black hole mass distributions based on the 
Press - Schechter theory of the structure formation in the Universe with quasar luminosity func- 
tions resulting from a physics-based emission model that takes into account the time-dependent 
phenomena occurring in the accretion disks. Quasar evolution and CDM models are mutually 
constraining, therefore our model gives an estimation of the exponent, n, of the power spectrum, 
P(fc), which is found to be —1.8 < n < —1.6. We were able to reject a generally assumed hypoth- 
esis of a constant ratio between Dark Matter Halo and the Black Hole mass, since the observed 
data could not be fitted under this assumption. Wc found that the relation between the Dark 
Matter Halos and Black Hole masses is better described by Mbh = ^dmh- This model provides 
a reasonable fit to the observed quasar luminosity function at redshifts higher than ~ 2.0. We 
suggest that the disagreement at lower redshift is due to mergers. Based on the agreement at 
high redshift, we estimated the merger rate at lower redshift, and argue that this rate should 
depend on the redshift, like (1 -I- z)^. 



Subject headings: accretion: accretion disks - cosinology:theory - quasars: general 



1. Introduction 

Since their discovery by Schmidt in the early 
60's, quasars have been an intriguing phenomenon 
and the mystery surrounding them is still far 
from being resolved. Soon after being identified 
as galactic nuclei, accretion onto a supermassive 
black hole became the leading model for power- 
ing these objects (Lynden Bell, 1969), whose lu- 
minosities can far exceed those of galaxies, but 
whose space density is around two orders of mag- 
nitude lower than that of bright galaxies. The 
question of whether quasar activity was a short- 
lived phenomenon that occurred repeatedly during 
a galaxy's lifetime (Cavaliere & Padovani, 1989; 
Siemiginowska & Elvis, 1997), or whether quasars 
never "wake up" again after a single, longer-term. 
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activity cycle (Haiman & Loeb, 1998; Haehnelt, 
Natarajan & Rees, 1998), leaving behind massive 
black hole remnants, has not yet been resolved. 

The quasar population clearly evolves with cos- 
mic time. The quasar space density above a 
given luminosity increases by a factor of ^ 50 
from the present to z ~ 2.5, then stalls and 
probably decreases to z '^ 5 (Hook, Shaver & 
McMahon, 1998). For many years the observed 
quasar evolution has been addressed phenomeno- 
logically, mostly via Pure Luminosity Evolution 
(PLE - Mathez, 1976) or Pure Density Evolution 
(PDE - Schmidt, 1968). Both descriptions have 
just one free adjustable parameter, but are phys- 
ically rather abstract. Moreover PDE and PLE 
are observationally indistinguishable when only a 
power law is observed, as is common. Instead 
physically based models may be constrained by 
a simple power law, since neither the slope nor 
the normalization is now arbitrary. An early at- 
tempt to explain the observed quasar space den- 



sity using a physics-based model was made by Ef- 
stathiou and Rees (1988), who used the Press - 
Schechter formahsm (Press & Schechter, 1974) for 
dark matter halo formation in order to predict the 
space density of the more luminous quasars. This 
work has been followed by several such attempts, 
with results that do not always agree (Cavaliere & 
Padovani, 1988; Haehnelt & Rees, 1993; Yi, 1996; 
Sicmiginowska & Elvis, 1997; Haiman & Menou, 
2000; Haehnelt, Natarajan & Rees, 1998; Salucci 
et al, 1999). Cavaliere & Vittorini (1998) sug- 
gest that the steep rise of the number density of 
quasars from early epochs to a redshift of ^^2.5 is 
due to the newly formed galactic halos (density 
evolution) while their fall at lower redshifts is due 
to merging processes (luminosity evolution), cre- 
ating a link between the phenomenological models 
and the more robust theory. Quasars are now a 
common ingredient of N-body simulations, semi- 
analytical or analytical models, which follow step 
by step the formation and evolution of galaxies 
within Dark Matter Halos and the merging pro- 
cesses (e.g. Richstone et al., 1998; Kauffmann & 
Haehnelt, 2000; Monaco at al., 2000). All of the 
above suggests that galaxy and quasar, both for- 
mation and evolution are closely related, not least 
because there is so far no evidence of quasar ac- 
tivity outside galaxies. 

Previous work derives the black hole mass dis- 
tribution from Cold Dark Matter (CDM) mod- 
els (Haehnelt & Rees, 1993; Haehnelt, Natara- 
jan & Rees, 1998) or from observations of the lo- 
cal non-active galaxies (Salucci et al., 1999), but 
then considers an essentially arbitrary exponen- 
tial quasar light curve: a short active phase with 
luminosities close to the Eddington limit at the 
beginning of the black hole's life, followed by a 
rapid exponential decay, due to fuel exhaustion 
caused by the high accretion rates needed. Such 
a quasar light curve is only an approximation of 
the accretion process and does not include the 
physics of the accretion disk which is thought 
to be formed during accretion onto a supermas- 
sive black hole in quasars (e.g. Lyndcn-Bell 1969, 
Shields 1978, Rees 1984). Accretion disks have 
been studied in detail in many types of binary sys- 
tems, where they usually show complex variability. 
Both observations and theory point to strong sim- 
ilarities in the behavior of accretion disks around 
Galactic black hole X-ray sources and of those 



around quasars (Tanaka & Lewin, 1995; Fiore 
& Elvis, 1997; Czerny, Schwarzenberg-Czerny & 
Loska, 1999). These similarities suggest quite a 
different quasar light curve. Large amplitude in- 
stabilities typically appear in accretion disks, lead- 
ing to repetitive outbursts as large as a factor of 
10^ on the timescale of thousand to million years 
depending on the disk model (Siemiginowska, Cz- 
erny and Kostyunin, 1996, Lin & Shields 1986). 
In the case of quasar black holes, these variations 
alone inevitably create a broad luminosity func- 
tion (Siemiginowska & Elvis, 1997). 

Unstable accretion leading to episodic out- 
bursts has a valuable feature. In our previous 
paper (Siemiginowska & Elvis, 1997) quasar light 
curves were derived by assuming a continuous 
fuel supply at a low constant accretion rate of 
~ O.OAMEdd ( where Msdci is the critical accre- 
tion rate corresponding to the Eddington luminos- 
ity for a given mass). At this rate a 10^ Mq black 
hole accretes some 10~'^MQ/yr while a 1O^M0 
black hole accretes IMq/yt, so a black hole would 
need 10^" years to change its mass by an order of 
magnitude. Consequently the problem of creat- 
ing hypermassive {IO^^Mq Soltan 1982, Phinney 
1983) quasar remnants does not occur, if the high- 
est initial black hole mass is 1O^M0. 

In this paper we build on the Sicmiginowska 
& Elvis (1997) work by combining the CDM de- 
rived black hole mass spectrum with the luminos- 
ity function for a particular mass derived from the 
time-dependent accretion disk theory. This results 
in an end-to-end model for quasar evolution that 
is completely physics based, excepting only the ra- 
tio of black hole mass to dark matter halo mass. 
Since this model links CDM and quasar evolution, 
the parameters of CDM constrain quasar evolu- 
tion, and the observed quasar evolution constrains 
CDM parameters. We examine some of the impli- 
cations of this model. 

2. Model and assumptions 

2.1. Accretion onto supermassive black 
holes 

We assume that each quasar hosts a super- 
massive black hole surrounded by an accretion 
disk. The disk is the main emission component 
and driver of the quasar activity. We calculate 
a quasar light curve assuming that quasar cmis- 



sion is due to a non-stationary accretion disk, thus 
the quasar luminosity during evolution of accre- 
tion disk is calculated by integrating the disk emis- 
sion at each epoch. Our time-dependent accretion 
disk model (Siemiginowska et al., 1996) takes into 
account thermal-viscous instabilities due to Hy- 
drogen ionization, which can develop in the disk 
(Smak, 1982; Lin & Shields, 1986). These insta- 
bilities are observed in both cataclysmic variables 
and Galactic Black Hole Candidates (Smak, 1984; 
Tuchman, Mineshige & Shields, 1990), where they 
cause factor of 10^ luminosity outbursts. Similar 
amplitudes are expected in the disks around super- 
massive black holes but the timescales are of or- 
der of hundred to hundred million years (Burderi 
et al 1998, Siemiginowska, Czerny & Kostyunin, 
1996; Mineshige & Shields 1990). As a conse- 
quence a single source can emit within a wide lu- 
minosity range depending on the current state of 
the disk. Even if there is a continuous fuel sup- 
ply the source will exhibit many active and qui- 
escent states, making 'quasar activity' a recurrent 
phenomenon occurring numerous times during a 
galaxy's lifetime. For example, for a black hole 
of IO^Mq, 5 outbursts occur in ^ 3.2 x 10® years 
(figure 1). Siemiginowska & Elvis (1997) derived 
the quasar luminosity function by integrating the 
time spent by a quasar in each luminosity state 
and assuming that this distribution is applied over 
a whole (uncoordinated) population with a range 
of black hole masses. 

We have checked a simplifying assumption in 
Siemiginowska & Elvis (1997) that the distribu- 
tion of time in each state is constant for all black 
holes masses. In figure 1 we show the light curves, 
due to emission from the accretion disk which ex- 
hibits thermal- viscous instabilities, for black holes 
with masses 10*^ Mq, IO'^Mq, IO^Mq and IO^Mq. 
The overall form of the light curves is the same, 
but for larger masses, outbursts are rarer than for 
lower masses, the high luminosity phases last much 
longer and the mean luminosity increases too. Fig- 
ure 2 quantifies this by showing how the number 
of outbursts per 10® yr and the duration of each 
outburst vary with the black hole mass. Both the 
number of outbursts and their duration depend 
quadratically on the black hole mass. 

Note that at 10® M© an outburst occurs every 
century, which means that in a sample of 100 AGN 
we should be able to observe one transition per 



year. The Piccinotti et al (1982) 2-10 keV se- 
lected sample has ~ 20 non-blazar AGNs, roughly 
equally divided between broad and narrow- line ob- 
jects. In the ~ 20 years that these sources have 
been studied, one (NGC 7582) has changed from 
purely narrow to broad-line status, followed by 
a considerable change in its luminosity (Aretx- 
aga et al., 1999). This is consistent with the ex- 
pected statistics. Monitoring of large (N~ 1000) 
emission-line selected AGN samples for a few years 
would strongly test this key premise of the model. 
Integration of the light curves over time shows 
that each source spends the greater part of its life- 
time (around ~ 70%) in a "quiescent" state and 
~ 30% in an active state (L > O.QOlLEdd) with ^ 
10% in a very active state {L > O.lLEdd)- These 
percentages are almost identical for all black hole 
masses, as shown in table 1. Deviations are no 
larger than 5% over a range of 4 orders of magni- 
tude, which is insignificant. The fraction of time 
a source spends in each activity state is essentially 
indeed independent of its mass, as assumed by 
Siemiginowska & Elvis (1997). 

2.2. Dark matter halos and black holes 

To derive the black hole mass function, we use 
the Press - Schechter (1974) formalism, a sim- 
ple yet quite adequate approach, which has be- 
come the basis of most formation scenarios of 
dark matter halos in CDM models (e.g. White & 
Frenk, 1991; Haehnelt & Rees, 1993). The Press 
- Schechter formalism predicts the mass fraction 
of the Universe contained in virialized structures, 
resulting from a gravitational hierarchical process 
(dark matter halos). The mass function of the 
dark matter halos (DMHs) can be written as: 



n{M, z)dM = - 



2 \ p 5c da 
M'^dM 



exp 



4)«M 

2a^J 

(1) 
where p is the comoving density. This formula 
gives the number of dark matter halos of mass M 
formed up to a redshift z. 5c, the critical over- 
density in linear theory for spherical collapse, has 
been considered constant, equal to 1.69 (e.g. Sahni 
& Coles, 1995), although it depends weakly on 
f^o a-nd A. The linear theory rms mass density 
fluctuation in spheres of mass M at redshift z is 
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Fig. 1. — Quasar light curves (accretion disk luminosity vs time in 10^ yr) due to thermal- viscous instabilities 
in the accretion disk around a central black hole. From top to bottom and left to right the masses of the 
black holes are IO^Mq, lO^M©, IO^Mq and IQ^Mq. Accretion rate at the outer edge of the disk is equal to 
O.QAMEdd of critical (Eddington limit) accretion rate. 



where 

D,^{1 + 3/2^o)(l + 3/20o + 5/2noz)-^ 

is a simplified form of the linear growth factor 
(Peebles, 1980) that represents the growth of a 
perturbation in the linear regime. For the calcula- 
tion of (T we used the approximation of the present 
rms mass density fluctuation in spheres of mass M 
(e.g. Fanet al., 1997): 



(To ~ CTS 



Ms 
M 



where Mg — 1.19 x 10 JIq is the mean mass in a 
sphere of radius Sh^^Mpc, a = (n + 3)/6 where n 



is the exponent of the power spectrum P{k) ex /c", 
and the corresponding erg = 0.67. Comparing to 
other results (Haehnelt, Natarajan & Rees, 1998) 
we verified that the shape of the mass function for 
masses greater than IO^Mq is not affected by this 
approximation. 

The DMH mass function for the standard CDM 
model, for n = — 1 and for various redshifts be- 
tween 5 and is illustrated in figure 3a. A single 
cosmological model (qo=0.5) will be used hence- 
forth. 

Equation 1 refers to halos harboring both el- 
liptical and spiral galaxies. Quasars are found in 
both elliptical and spiral galaxies (McLeod et al., 
1995; McLeod & Katris, 1995; Taylor et al., 1996; 




Fig. 2. — (a) Variation of the number of outbursts (high luminosity active state) per 10^ yr and (b) duration 
of the outbursts (high luminosity state) (in 10^ yr) versus the black hole mass. The solid lines illustrate the 
linear fits. 

Table 1: Fraction of time sources of different mass spent in the different activity states 
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Bahcall et al., 1997; Boyce et al., 1998). Narrow 
lined (type 2) AGN can also reside in both early 
type galaxies and in spirals (e.g. Salucci et al., 
1999). We assume that the distinction between 
quasars and Seyferts is purely one of luminosity 
(we also include both radio-loud and radio-quiet 
AGN in term 'quasar'), the latter occurring in spi- 
ral galaxies, since we consider that they are both 
the manifestation of the same basic phenomenon. 
The black hole mass function can be directly de- 
rived from the dark matter halo mass function. 
We assume that in each newly formed dark matter 
halo a black hole is "instantly" created, i.e. that 
all galaxies harbor a black hole. This hypothesis 
has been supported by several observations (e.g. 
Keel, 1983; Ho et al., 1997; Magorrian et al., 1998; 
Richstone et al., 1998; Ford et al., 1998). Haiman 



& Loeb (1998) assumed that black hole formation 
is restricted to halos whose masses are bigger than 
1O*((1 + z)/1O)-3/2M0 with the mass of the black 
hole being a constant fraction of the dark matter 
halo mass: Mbh = W^Mumh, where a =-3.2 
for all black hole masses. We adopt this mass re- 
lation as a baseline, but examine the possibilities 
of varying a between -4.5 and -3.2. The result- 
ing black hole mass functions are shown in figure 
3b. Values close to -3.2 are consistent with the 
MsH/Mbuige = 0.006 of Magorrian et al. (1998), 
while values close to -4.5 are consistent with the 
MsH/Mbuige = 10-3-5, given by Wandel (1999), 
where M^uige denotes the mass of the galactic 
bulge. We also consider the possibility of an a that 
is a function of the DMH mass, (3 x log{ME>MH), 
which is equivalent to a non-linear relation be- 



tween the BH and DMH masses, Mbh — ^dmh- 
The dash-dotted hne in figure 3b traces the black 
hole mass function for /3 = —0.332, a value cho- 
sen by fitting the observations, as explained in the 
next section. 

3. Deriving the luminosity function 

The luminosity function is the product of the 
fraction of time each source of a given mass spends 
in each luminosity (activity state) and their space 
density (Siemiginowska & Elvis, 1997). More for- 
mally, the luminosity function of a population of 
quasars is the result of the convolution of their 
mass function with the light curve, L(t), of an in- 
dividual source of each mass: 



$(i,z)= <^{L,z,M)n{z,M)dM, 



(2) 



where n{z, M) is the mass function and $(i, z, M) 
is the fraction of objects that contribute to a given 
luminosity at a given redshift (equivalent to fon in 
equation 3). 

This is equivalent to the Haiman & Menou 
(2000) suggestion that the observed quasar lumi- 
nosity function dNobs{z, M)/dL, can be directly 
compared with the Press - Schechter DMH mass 
function, dn{z,M)/dM, to derive the luminosity 
of each halo of mass M: 



d<i>obs{z,MBH) dL 



dL 



dM 



^.fo 



dn{z,M) 
' dM '' 



(3) 



where fon is the quasar "duty cycle" defined as the 
fraction of black holes that are active at a given 
redshift, z. Note that equation 2 is the integrated 
form of equation 3, if the halo mass function is 
replaced by the black hole mass function. 

By integrating the light curves over timescales 
that allow many outbursts to occur, we calculated 
the time a source spends in the different activity 
states for a range of black hole masses. This is 
equivalent to calculating, for a certain redshift, 
z, the fraction of sources of a given mass, M, 
in a given activity state, L, <I>(L,z,M). All the 
above are illustrated in figure 4a, where we show 
the contribution of sources with black hole masses 
from 10® Mq to lO^M©, to the overall luminos- 
ity function, for a Press - Schechter derived black 



hole mass function at a redshift of 2. The sources 
with black hole masses within the abovementioned 
range have luminosities within the observed lim- 
its (10'*^erg/sec < L < W^^erg/sec). Masses 
have been converted to luminosities through the 
light curves and a was held constant, Mbh — 
10~^-^M]jMH, for all black hole masses. 

Note that in each case, a low luminosity state 
for a high mass source, is a high luminosity state 
for a lower mass source, and vice versa. Fig- 
ure 4b shows that at each redshift, the slope de- 
fined by the high activity sources gives roughly 
the shape and slope of the luminosity function. 
In the case of an upper or lower mass limit, the 
luminosity function will follow the shape of the 
corresponding light curves, creating two breaks: 
one towards the lower and one towards the higher 
luminosities. Figure 5 illustrates the luminosity 
function of the form ^{MB)dMB, for different val- 
ues of the parameter n (Fig. 5a) and the pa- 
rameter a (Fig. 5b), including black hole masses 
IO'^Mq < Mbh < IO^Mq, after converting the 
luminosities into absolute magnitudes, Mb- The 
shape of the luminosity function and its evolution 
depend on these two parameters. Our discussion is 
focussed on the effect of varying these parameters. 

3.1. Effects of changing CDM Power Spec- 
trum Slope (n) 

Here we examine the effect of changing the val- 
ues of the exponent n of the CDM power spec- 
trum. For a CDM- type spectrum P{k) ex fc", n 
takes values in the interval —2 < n < — 1 on scales 
of ~ 8h~^ Mpc (e.g. Fan et al., 1997). For many 
years, the galaxy correlation function data were 
consistent with n = — 1.2 (Peebles, 1980; Davis & 
Peebles, 1983). However, new surveys show that 
the shape of the power spectrum has a more com- 
plicated form and/or that the values of n vary be- 
tween small and large scales (e.g. Lin et al., 1996; 
Hermit et al., 1996; Le Fevre et al., 1996; Tadros et 
al., 1999). We do not attempt to introduce a more 
complex form for the power spectrum. However, 
we do find that varying the exponent, n produces 
strong effects on the quasar luminosity function. 
Observations show that the quasar space density 
increases as $(z) = <l>o(l-|-z)'^ from the current 
epoch up to a turnover redshift, zt, of 2-3, where 
it stops or possibly drops (e.g. Hook et al., 1998). 
Our models naturally produce a zt by integrating 



■^ 


2 


,:^. 




(q) n = -l - 


--^ 





^"^^^""^^ 




- 


^ 




"~~^^ 


%< 


' 












(J' 







"^■- 


~^r*%^ 


'"^ 


-2 






"~"^^^. - 











'\\-- 


Q_ 








■ \^''^- 


5 








\ \\ ■^•." 


^-^ 


-4 


" z=5 




■■ \ \ - 


I 




■ 2=3 




\ \" 


2 








\ - 


^ 


-6 


z = z 




\~ 


Oi 




■ -.-. z=0 




■, \- 


- 


p 






\ - 



10 11 12 

log (M„.„) 



14 



01 




-6 



2 

19^ 



Qi 




(b) z=3 



a = -3.2 
a = -5.8 
a = -4.5 
a = -0.332*log( 



log (Mbj 



Fig. 3. — (a) Press - Schechter DMH mass function for z=5 (dotted line), 3 (solid line), 2 (dashed line) and 
(dash-dotted line) for the cosmological model qo=0.5 and for n ~ —1. (b) Black hole mass functions for 
z = 3. The solid line corresponds to a~-3.2, the dotted line to a^-3.8, the dashed line to a=-4.5 and the 
dash-dotted line to a = —0.332 x log{MoMH)- 



over the whole redshift range the predicted lumi- 
nosity function, for sources with absolute magni- 
tudes brighter than Mb = — 26. We find that the 
Zt predicted by our models is a parabolic function 
of n. Figure 6a shows the logarithmic normalized 
space density of quasars brighter than Mb = — 26 
as a function of z, for 4 different values of n (with 
a — —4.5). The peak of the distribution, zj, 
changes with n: the steeper the power spectrum 
the lower the zt- Figure 6b shows the relation be- 
tween Zt and n, which is roughly parabolic. More- 
over, this relation is almost independent of the 
value of a making this a robust result, since a 
is the only physically undetermined parameter in 
the model. 

The observed Zt is at least as high as 2.5 (Hook 
et al., 1998) and may be as high as ~ 4 (Miyaji 
et al., 2000), but quasar space density declines 
rapidly for 2; > 4, as suggested by both optical and 
radio observations (Kennefick et al., 1997; Hook et 
al., 1998) and theoretical models (e.g. Cavaliere & 
Vittorini, 1998, give a zt around 3.5). This puts 
a strong constraint on n, which can only take the 
values between -1.6 and -1.8 (figure 6b). These 
values are in excellent agreement with the values 



given by Blanchard et al., 1999, derived from the 
galaxy clusters. We restrict the rest of our study 
to this range of n. 

3.2. Dependence of Halo to black hole 
mass ratio (a) on mass. 

In the next two sections we discuss how the pa- 
rameter a affects the luminosity function and we 
consider two possibilities: (1) a is constant for all 
masses; (2) a is related to the mass of the DMH. 

3.2.1. Constant (a) 

Most theoretical models have assumed for sim- 
plicity that the black hole mass is a constant frac- 
tion of the dark matter halo mass, independent of 
the masses involved (e.g. Haiman & Loeb, 1998; 
Haiman & Menou, 2000). Observationally the re- 
sults of Magorrian et al (1998) suggest that this is 
roughly obeyed. Therefore we first examine this 
possibility and compare our predictions with ob- 
servations, both for low and for high redshifts. 
Figure 7 illustrates the predicted luminosity func- 
tion in the case of a = —3.2 and a = —4.5 for 
different values of n ranging from -1.6 to -1.8, at 







,■& -7 







^ 


, 1 , , , 1 , , , 1 , , , 1 , , , 1 , _ 




: \ (q) : 


^ 


: v^.,_i_^ ^ 


t; 


: ^r^,..,...„ \ 


6 






1 .-'-'''' ^ : 


/ 


\ : 


Q 





40 42 44 46 48 

bg(L) (erg/sec) 







,^ -7 







-8L 




40 42 44 46 

bg(L) (erg/sec) 



48 



Fig. 4. — (a) Contribution of the different black hole mass sources to the luminosity function at z — 2, for 
black hole masses ranging from IO^Mq (solid line) to 10^ Mq (dash-dotted), a was constant and equal to 
-3.2 for all black hole masses, (b) The slope defined by the sources at high luminosity state gives roughly 
the slope of the luminosity function. In the case of an upper (here 10^ Mq) or lower (here IO^Mq) mass 
limit, the luminosity function will follow the shape of the corresponding light curves, creating two breaks. 



z = 3.25. The symbols are the observed luminos- 
ity function from Warren et al. (1994). 

In both cases, the slope of the predicted lu- 
minosity function is different and much smoother 
than the observed luminosity function and no 
break is present. The observed luminosity func- 
tion cannot be fitted with the constant a model. 
For a high value of a (-3.2) the luminosity func- 
tion is over-estimated at high luminosities, while 
for a low value {a = —4.5) the numbers of faint 
quasars are under-predicted. The same tenden- 
cies are seen at all redshifts. This suggests that a 
cannot be constant, but may be a function of the 
masses involved. 

3.2.2. a dependent on the Mdmh 

There is no a priori reason that the black 
hole mass and the dark matter halo mass should 
have constant ratio, so in this section we dis- 
cuss the effects of varying a. We assume that 
a depends on the DMH's mass with the relation 
a = (3 X log{MoMH), or, in other words, that 
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In order to obtain the best 



value of a parameter (3 we fit this function to the 
observed luminosity function from Warren et al. 
(1994) at redshifts 2.6 and 3.25, shown in figure 8. 
We found that (i = -0.332to;ot6 (or 1+/? = 0.668) 
gives the best fit, which is the value used below. 

The black hole masses are now defined from the 
relation: Mbh = ^dmh- This relation gives a=- 
3 for Mbh = IO^Mq and a=-4.4 for Mbh = 
10^ Mq, which are in agreement with the observed 
values (Magorrian et al., 1998 and Wandel, 1999). 

Figure 8 illustrates the predicted luminosity 
function for z = 2.6 and z = 3.25 for 3 values of n 
from -1.6 to -1.8. The symbols show the Warren 
et al. (1994) data. The fit is much better than 
for any constant value of a. The bright end of our 
predicted luminosity functions at z = 2.6 seems to 
be over-predicted, when compared to observations 
(figure 8a) . Alternatively, the observed bright end 
may have been under-estimated, as suggested by 
recent surveys, in particular the Hamburg/ESO 
survey (Wisotzki et al., 1996, Kochlcr et al, 1997 
and Wisotzki, 2000). We conclude that a must 
be a function of the DMH in order to properly 
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Fig. 5. — Predicted luminosity funetion for z 
n is constant. 



2. (a) n varies while a is kept constant, (b) a varies while 



describe the most recent observations. 

Our physical model leads to changes in the 
luminosity function equivalent to either PLE or 
PDE, up to a redshift of ^ 4, because its form is 
close to being a simple power law. 

3.3. Comparison with observations - dis- 
cussion 

The predicted evolution fits the data quite well 
at high redshifts, e.g z > 2.0. For redshifts below 
^ 2, however, the predicted luminosity function, 
in terms of PLE with L ~ L={1 + z)^p evolves as 
7p = 0.5. This disagrees with the observed lumi- 
nosity function issued from the following samples: 
AAT (Boyle et al., 1991), HES (Wisotzki et al., 
2000), LBQS (Hewett et al., 1995) and the CFRS 
quasars (Le Fevre et al., 1995), for redshifts 0.5 
to 2.0 (figure 9), that evolves much more strongly 
with, 7o ~ 3.5, although the slopes of the two 
(observed and theoretical) luminosity finctions are 
the same, as seen in figure 9. In this figure, the 
symbols show the observed luminosity function 
issued from the abovementionned samples while 
the dashed, dotted and thick solid lines illustrate 
the predicted luminosity function for n = —1.8, 
n = — 1.7 and n = —1.6, respectively. 



For many years, the best established form of the 
observed quasar luminosity function was a bro- 
ken power law, first given for optically selected 
quasars with redshifts 0.3 < z < 2.2 by Boyle 
et al. (1988), confirmed by La Franca & Cris- 
tiani (1997) for other optical surveys and by Page 
et al. (1996) for X-ray selected AGN. However, 
recent surveys, in several wavelengths, show that 
this assumption was probably due to selection ef- 
fects. For z > 2.2 (as Lya shifts from the U to 
the B-band) multicolor selection techniques result 
in biased samples and favor the redshift regions 
at which quasar colors make them distinguishable 
from stars. Giallongo and Vagnetti (1992) argued 
that the observed break was due to the uncertain- 
ties on the K-correction due to the dispersion in 
quasar optical slopes. This seems to be confirmed 
by ROSAT AGN data (Miyaji et al., 2000). More 
recent surveys, which account for this effect, such 
as the Hamburg/ESO Quasar Survey (which uses 
emission lines in addition to colors, Wisotzki et 
al, 1996, Koehler et al, 1997 and Wisotzki, 2000) 
show that the number of bright objects was under- 
estimated and that the combined local luminosity 
function for Seyferts and quasars could be fitted 
with a single unbroken power law, as argued by 
Giallongo and Vagnetti (1992). Pei (1995) also 
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suggests that a single power law of the form L^/^ 
could fit the data as well as a broken power law. 
This single power law is represented by the thin 
solid line in figure 9. In terms of luminosity, galax- 
ies with central black hole masses of the order of 
l(fMQ will be seen as Seyfert galaxies, while the 
higher masses will result in quasars (figure 4). 

Even if the observations are altered somewhat 
additional physical effects are likely to be impor- 
tant, since we derived the mass function using only 
the Press - Schechtcr DMH mass function. For ex- 
ample galaxy mergers could significantly modify 
the assumed mass function (e.g. Haiman & Loeb, 
1998; Kauffmann & Hachneh, 2000). Mergers are 
rare in the early Universe but start becoming im- 
portant at a redshift around 2, an epoch of massive 
galaxy formation (Kats et al., 1992; Kontorovich 
et al., 1992). Only mergers that result in the for- 
mation of a single object will cause a decrease of 
the space density of the black holes. The black 
holes, due to the dynamical friction, will soon be 
found in the center of the newly formed object 
and merge (e.g. Binney & Tremaine, 1994), thus 
influencing the quasar space density. If this is the 



primary cause of the low z disagreement of our 
model with observation then the difference of the 
observed to predicted evolution rates, 7o — 7p, pro- 
vides a crude estimate of the merger rate, 7™, for 
redshifts between 0.5 and 2, i.e. 7™ — 3. Hence 
the merger rate can be expressed as a function of 
the redshift: ~ (1-1-2;)^. Interestingly this value is 
in agreement with several other results based on 
observations: Zepf & Koo (1989) studied a com- 
plete sample of 20 close pairs of faint (B < 22) 
galaxies and found that the frequency of interact- 
ing galaxies increase as {l + zY'^^"^'^ . Burky et al. 
(1994) used four fields observed with the Hubble 
Space Telescope and found that the galaxy merger 
rate increases with redshift as (1 -I- z)^-^^^-^ . Pat- 
ton et al. (1997) used a sample of 545 field galaxies 
and estimated that the merger rate changes with 
redshift as (1 + z)2-8±o.9_ 

Note that less extreme galaxy interactions, 
where the objects remain independent entities 
afterwards will not influence our results, if we 
assume that whatever happens in the vicinity of 
an object will not influence its accretion rate and 
therefore its light curve. In reality, such interac- 
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Fig. 8. — Predicted luminosity function for z = 2.6 and 3.25, for n—-l.6 (solid line), -1.7 (dotted line) and 
-1.8 (dashed line). Symbols show the observed luminosity function of Warren et al. (1994). 



tions are likely to affect the accretion rate and will 
be taken into account in the future work. 

The disagreement between the observed and the 
predicted luminosity function at z < 2 may also 
be due to the role of obscured (type II) AGN. The 
estimates of their number density are not yet accu- 
rate. Recent work (Salucci et al., 1999) estimates 
the local mass function of the dormant black holes 
within elliptical galaxies, assuming that the nu- 
clear activity is a single short event. Salucci et 
al., however, do not take merger effects into ac- 
count, introducing uncertainties in the calculated 
massive dark object mass function. 

We used a black hole mass distribution with 
black holes ranging from ^ IO^Mq to IO^Mq. Ac- 
cording to our model, accretion onto such black 
holes can produce luminosities within the observed 
range lO'^i <L< lO^^ erg/sec (10"'* < L/Lsdd < 
1). If we allow for a higher black hole mass cut-off 
at ~ 10^ Mq we could obtain the break in the lu- 
minosity function. However, there is evidence for 
black holes with lower masses (e.g. Reid et al., 
1999, Sag A* - 2.6 x IO^Mq). Therefore a break 
cannot be imposing a simple cut-off mass. 

Our model gives a consistent physical connec- 
tion between high and low luminosity sources and 



does not predict a break in the luminosity func- 
tion for the black hole mass distribution within 
the observed mass range. 

The model presented here can be extended and 
used for the study of the UV-background and 
the re-ionization history of the Universe (Siemigi- 
nowska et al., 2000, in preparation). Issues such 
as whether quasars are entirely responsible for the 
re-ionization of the Universe or whether Popu- 
lation III stars are required (see Madau 1999), 
the redshift at which it occurred, whether it is 
homogeneous or inhomogeneous or its effects on 
the primordial power spectrum, can be examined. 
The results are likely to be different from previous 
attempts because the emission spectrum derived 
from a non-stationary accretion disk is different 
from the one predicted by a stationary model at 
the UV wavelength range (Siemiginowska et al., 
1997). 

4. Conclusions 

In this paper we presented a theoretical model 
for calculating the quasar luminosity function. 
This model, unlike previous ones, takes into ac- 
count the physical processes occurring in the ac- 
cretion disk around a black hole. Black holes are 
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assumed to form instantly within all dark matter 
halos, with masses related to the halo mass as: 

We first restricted the range for the two free 
parameters of the model: the power spectrum 
exponent, n and the ratio of black hole mass 
to DMH mass, a. For any reasonable a value, 
(—4.5 < a < —3.2), n has to be within the nar- 
row range -1.6 to -1.8 in order to match the ob- 
served redshift of the peak of quasar evolution. 
We compared the predicted luminosity function to 
the data and found that a non-linear relation be- 
tween the BH and DMH masses {Mbh = Af^"^^) 
gives a much better approximation than a linear 
relation {Mbh = ^0"Mdmh)- Our model pre- 
dicts no break in the luminosity function, assum- 
ing quasars and Seyfert galaxies to be objects of 
the same kind but of different luminosities. Ob- 
servationally this model is then equivalent to both 
PLE and PDF, but physically has quite different 
implications. 

No assumptions have been made on the merger 
processes, although they impact the quasar evolu- 



tion. Mergers that result in the black hole fusion 
(and consequently in a decrease of the black hole 
space density) will also influence the luminosity 
function at redshifts lower than ^ 2, when mergers 
become important. We estimated that the merger 
rate should depend on redshift, (l-|-z)^, consistent 
with direct measurements. 

Built into our model is the idea that quasar-like 
activity is a recurrent phenomenon. Practically all 
galaxies go through that state many times during 
their lives. Recent observations support the idea 
that black holes reside in all galaxies. 

This work was supported by NASA grants 
NAG5-3391, NAG5-4808, NAG5-6078, and NAS8- 
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Fig. 9. — The predicted quasars luminosity function for -1.8 < n < -1.6 at (a) z — 0.5, (b) z — 1.0, (c) 
z = 1.5 and (d) z = 2.0 (thick hncs). The symbols represent the observed luminosity function for the samples 
AAT (Boyle et al., 1991), LBQS (Hewett et al., 1995), HES (Wisotzki et al., 2000) and CFRS (Lc Fcvrc et 
al., 1995) (Wisotzki, private communication). The thin solid line represents the single power law LF of Pei 
(1995). 



15 



